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1 Introduction 



The generalization of the concept of derivative and integral to a non-integer 
order a has been subjected to several approaches and some various alternative 
definitions of fractional derivatives appeared [1-6]. In the last few years 
fractional calculus was applied successfully in various areas , e.g. chemistry 
, biology, modelling and identification, electronics and wave propagation. 
Fractional calculus , has played an important role in engineering, science, 
and pure and applied mathematics [7-9]. Fractional derivatives were applied 
in recent studies of scaling phenomena [10-12]. Classical mechanics is one of 
the fields when fractional calculus found many applications [13-19]. Riewe 
has used the fractional calculus to obtain a formalism which can be applied 
for both conservative and non conservative systems [13, 14]. Although many 
laws of nature can be obtained using certain functionals and the theory of 
calculus of variations, not all laws can be obtained by this manner. As 
it is known, almost all systems contain internal damping, yet traditional 
energy based approach cannot be used to obtain equations describing the 
behavior of a non conservative system [13, 14]. Using the fractional calculus 
one can obtain the Lagrangian and the Hamiltonian equations of motion for 
the nonconservative systems . 

The understanding of constrained dynamics j20], both at the classical 
and quantum level, has been a subject of long standing theoretical interest, 
which has seen important contributions ever since Dirac's quantization of the 
electromagnetic field. The path integral approach and the canonical one are 
two main approaches of quantization. 

Recently, an extension of the simplest fractional problem and the frac- 
tional variational problem of Lagrange was obtained [17,18]. Even more 
recentlythis approach was extended to Lagrangians with linear in velocities 
[21,22] , which represents a typical example of second-class constrained sys- 
tems in Dirac's classification [20]. These Lagrangians are important because 
their Euler-Lagrangian equations become systems of first order differential 
equations in contrast with second order corresponding to the regular ones. 
In addition, these systems may possess gauge symmetries and gauge ambi- 
guities. 

From these reasons it is interesting to study the fractional Hamiltonian 
of constrained systems. 

The aim of this paper is to obtain the fractional Hamiltonian equations 
of motion for Lagrangians with linear velocities. 
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The plan of our paper is as follows: 

In Sect. 2 some basic tools of fractional derivatives as well as Riewe's ap- 
proach of the fractional Lagrangian and Hamiltonians are presented. In Sect. 
3 the Euler-Lagrange equations were obtained using the Agrawal's approach 
and the fractional formulation of systems with constraints is introduced. In 
Sect 4 the fractional Hamiltonian of the systems possessing linear velocities 
is analyzed. Conclusion is presented in section 5. 

2 Fractional Lagrangian and Hamiltonian for- 
mulations 

In this section we briefly present the definition of the left and right derivatives 
together with Riewe's formulation of Lagrangian and Hamiltonian dynamics. 
The left Riemann-Liouville fractional derivative is defined as 

v ' a 

and the right Riemann-Liouville fractional derivative has the form 

where the order a fulfills n — 1 < a < n and V represents the Euler's gamma 
function. If a is an integer, these derivatives are defined in the usual sense, 
i.e., 

Now we shall briefly review Riewe's formulation of fractional generaliza- 
tion of Lagrangian and Hamiltonian equations of motion [13,14]. The starting 
point is the action function of the form 

S= f L({q r n ,Q r nf },t)dt. (4) 

J a 

Here the generalized coordinates are defined as 

q r n := {aD?) n x r (t), Q r n , := x r (t), (5) 
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and r = 1,2,..., it! represents the number of fundamental coordinates, n = 
0, N, the sequential order of the derivatives defining the generalized coordi- 
nates q, and n' — 1, N' the sequential order of the derivatives in definition 
of the coordinates Q. A necessary condition for S to posses an extremum for 
given functions x r (t) is that x r (t) fulfill the Euler-Lagrange equations [13,14] 

f)T N FIT N ' f)T 

ag n= i cg n n / =1 a(^ n / 

Using the references [13,14], the generalized momenta have the following 
form 



N 



E (W""- 1 ^ 



Pi 

fc=n+l 

Ar ' FIT 

<>= E ( a DD w -^. (7) 

k=n'+l 

Thus, the canonical Hamiltonian is given by 

R N-l R N'-l 

h = E E « + i + E E - (8) 

r=l n=0 r=l n'=0 

The Hamilton's equations of motion are as follows 

^ = -^- = (9) 
For n = 1, N, n' — 1, N' we have the following equations of motion 

^ 0L D^ + J)^. (11) 



The remaining equations are given by 



0H 4i = «, |^ = Q; +1 = ^q;„ (12) 



dp r n ^ IW a</ 

a// _aL 

where, n = 0, JV, n' — 1, AT'. 



(13) 
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3 Fractional Euler- Lagrange equations 

Recently Agrawal , have obtained the Euler-Lagrange equations for fractional 
variational problems [17]. In the following we like to present briefly his 
approach. 

Consider the action function 

%o\ ••••%*]= f L({q r n ,Q r n ,},t)dt, (14) 

J a 

subject to the independent constraints 

$ m (t, ql q r n , Q r n ,) = 0, m < R, (15) 
where the generalized coordinates are defined as 

q r n := ( a Bf) n x r (t), Q r n , := ( t B P b ) n 'x r (t). (16) 

Then, the necessary condition for the curves q^, q$ with the boundary 
conditions q^a) = q™ , q^b) = , r = 1, 2, , R, to be an extremal of the 
functional given by equation (14) is that the functions ql satisfy the following 
Euler-Lagrange equations [17]: 

f)T N r)T N ' r)T 

g-r+T, (Wg- + E (« D °)"'aoT = 0, (17) 

oq n= i oq n n t = i cy n / 

where L has the form [17] 
L({q r n , Q r n ,}, t, X m (t)) = L({q r n , Q r n ,}, t) + X m m m (t, ql..., q*, q r n , Q r n ,). (18) 
Here the multiple X m (t) G R m are continuous on [a, b]. 

3.1 Fractional Hamiltonian formulation of Agrawal's 
approach 

In order to obtain the Hamilton's equations for the the fractional variational 
problems, we re-define the left and the right canonical momenta as : 

N 8T 

k=n+l 0( ik 

N ' Fit 
<>= E (.D^-"'- 1 ^. (19) 

k=n'+l 
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Using (|19|).the canonical Hamiltonian becomes 

R N-l R N'-l 

^ = EE PWn + i + E E K'QWi - l- (20) 

r=l n=0 r=l n'=0 

Then, the modified canonical equations of motion are obtained as 

{q r n ,H} = t -D a bP r n , = a D t Q <„ (21) 

W ,H} = tV a b Po + aV> r o, (22) 

where, n = 1, N, n! = 1, N'. 

The other set of equations of motion are given by 



{Pn, H} = q r n+1 = aD?ql {<„ 8} = Q r n+1 = t D«Q r n „ (23) 
dH dL 
~dt ~ ~~dt' 



(24) 



Here, n = 0,...,N, n' = 1,...,N' and the commutator {,} is the Poisson's 
bracket defined as 

OAdB dB dA dA dB dB dA 

{ , " " + " (25) 

where, n = 0, N, n' = 1, N'. 

4 Equivalence of fractional Hamiltonian and 
Lagrangian formulations for systems with 
linear velocities 

Recently, for < a < 1, the Lagrangians with linear velocities were inves- 
tigated in [21]. For example, the Euler-Lagrange equations of the following 
Lagrangian 

L' = cij (V) aDtl j ~ V (V) , (26) 

were obtained as [21] 

—^- a T> t q J + t T> b a k (q) - — ^— = . (27) 
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Now we would like to obtain the Hamiltonian equations of motion for the 
same model. Let us define 



( a Bt) n q v n = 0,l,...,iV-l, j = 1, 2, R. (28) 



The generalized momenta are given by, 

The canonical Hamiltonian reads as 

H=(p° j -a j (x°))x) + V(x°). (30) 
The Hamiltonian equations of motion are calculated as 

^T=p5-^) = (*D?)^ = 0. (31) 

In fact, this equation is the primary constraint in the Dirac's formalism. The 
other equations of motion are calculated as 

dH_ da^) dV(x°) 

dxl~ dxl 3 d<l k ~ [tDb)Pk ' (32) 
f)H 

=x 1 k = ( a Pt)q k (33) 



dpi 

Making use of equation (31) we obtain 

-^- o D t ^+ t D b a fc ( g J--^- = 0. (34) 

It is obvious the equivalence between equations (34) and (27). An inter- 
esting point to be specified here is that, primary constraints in the Dirac's 
formalism are present as equations of motion in our treatment, while the other 
Hamiltonian equations of motion are equivalent to the Lagrangian equations 
of motion as given in [21]. 
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5 Conclusion 



One of the main problems encountered in applying the fractional calculus 
to a given singular Lagrangian and Hamiltonian is the existence of multiple 
choices of the possible fractional generalizations . In addition, the solutions 
of the fractional Euler-Lagrange equations contain more information than 
the classical ones. In this paper, Hamiltonian equations have been obtained 
for systems with linear velocities, in the same manner as those obtained by 
using the formulation of Euler-lagrange equations for variational problems 
introduced by one of us [21]. We study the general model for systems with 
linear velocities and it was observed that the Hamiltonian and Lagrangian 
equations which are obtained by the two methods are in exact agreement. 
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